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Abstract 

In aerospace computational fluid dynamics (CFD) 
calculations, the Delaunay triangulation of suitable 
quadrilateral meshes can lead to unsuitable triangu- 
lated meshes. In this paper, we present case studies 
which illustrate the limitations of using structured 
grid generation methods which produce points in a 
curvilinear coordinate system for subsequent trian- 
gulations for CFD applications. We discuss condi- 
tions under which meshes of quadrilateral elements 
may not produce a Delaunay triangulation suitable 
for CFD calculations, particularly with regard to high 
aspect ratio, skewed quadrilateral elements. 

1 Introduction 

In computational fluid dynamics (CFD) applications, 
the problem domain must be discretized into meshes 
(or grids) over which the governing equations of fluid 
dynamics are solved. The two major classes of grids 
for aerospace CFD applications are structured grids 
and unstructured grids. Structured grids sire curvi- 
linear grids designed so that the neighbors to any 
element are implicitly known. These grids have been 
studied for quite some time and techniques for their 
construction are well understood [2], Unstructured 
grids are composed of elements in which neighbors 
must be explicitly listed. The component elements 
are usually, but not necessarily, triangular. These 
grids are currently not as widely used in aerospace ap- 
plications, and have been the object of recent interest 
[3, 4, 5]. Several properties of Delaunay triangulation 
[6, 7, 8, 9] make it attractive to use in unstructured 
grid generation. However, a major drawback to this 
method is the need for a separate method of point 
generation; a straightforward approach to this draw- 


back is the use of a structured grid generator to create 
the necessary points [10, 11, 12, 13]. 

Certain features of structured grids are useful to 
maintain in unstructured grids. Most commonly, 
structured grids are body-fitted curvilinear meshes 
where contours follow the object boundaries. Struc- 
tured grids may also contain very high aspect ratio 
elements. This allows properties of the problems be- 
ing solved to be exploited; in general, very high so- 
lution gradients can exist perpendicular to surfaces, 
and very small solution gradients tangent to surfaces. 
Unstructured meshes should exhibit the same Struc- 
ture” which orients the cells along a feature of in- 
terest. When a structured grid contains high aspect 
ratio grid cells, the Delaunay triangulation of the grid 
points may not maintain the original grid lines of the 
structured grid (Figures 1 and 2). Orientation of the 
cells in a particular direction is lost when this occurs. 

This study investigates conditions under which De- 
launay triangulations of points from structured grids 
will maintain the original grid structure. The effects 
of skew and aspect ratio on the Delaunay triangula- 
tion Me studied. 

2 Preliminaries 

First consider a general case of points distributed ar- 
bitrarily along straight lines a distance h apart, with 
the points no more than s apart, and s > h. What 
relationship must exist between s and h so that the 
lines are guaranteed to be in the Delaunay triangula- 
tion of the points. 

Let Pi = (—s/2, 0) and P 2 = (s/ 2, 0) be two points 
on the line y = 0 (the results can be generalized to 
any location, but these points were chosen for simplic- 
ity of analysis). One scenario could have points on 
successive contours shifted by s/2, e.g., R x = (Q f h) 
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and Ft? = (0 ,—h). For the edge PiP 2 to be included 
in the triangulation of the points, the circumcircle for 
-Pi, P 2 , Ri must not contain P 2 , and likewise, the cir- 
cumcircle for Pi, P 2 , i ?2 must not contain R\. The 
circumcenter of Pi, P 2 , Pi is (0, h/2 — s 2 /8h). The 
circumcenter of Pi, P 2 , P 2 is (0,-/i/2 -f s 2 /8h). To 
guarantee that edge PiP 2 is in the triangulation, then 

-h/2 + s 2 /8h < h/2-s 2 /8h 

s 2 /h 2 < 4 

s/h < 2 

assuming distances positive. 

Without knowing anything else about the distribu- 
tion of points along the contour lines, we must con- 
clude that the aspect ratio of the triangles can be no 
greater than 2 to guarantee that all original edges are 
in the triangulation. The introduction of structure 
into the points locations allows for larger aspect ra- 
tios under certain conditions. The remainder of this 
section defines some of the concepts to be used, and 
following sections explore this idea in further detail. 

Definition 1 A structured grid is an undirected 
graph G=(V,E) with vertex set V and edge set E such 
that 

V = {v.-j I Vvij3(xij,yij) e £ 2 , 

0 < * < m, 0 < j < n} 

e = {(u t w) | vtnjev; 

1 ^ m => (vij , flj+ij) G E 

j n ^ ( v ij i v ij+i) € E } 

Definition 2 A graph G = (V t E) is called 
Delaunay-embeddable if G is a subgraph of the De- 
launay triangulation ofV, D(V). 

Definition 3 A rectangular structured grid is a 
structured grid G — (V,E) where the location of each 
Vij is given by: 

Xij = i • s 
Vij = j-h 

where s, h G 3?, 0 < i < m, and 0 < j <n. 

Points in a rectangular structured grid have a con- 
stant spacing in x and a (possibly different) constant 
spacing in y. 

Definition 4 The aspect ratio of a structured grid 
element is the ratio s/h , where s is the length of a 
long side and h is the separation between the sides of 
length s. 


Definition 5 A skewed structured grid with aspect 
ratio s/h is a structured grid where the location of 
each Vi t j is given by: 

Xi t j = * • s + j • dx 

Vij = i' h 

where dx < s/2, 0 < * < m, and 0 < j <n. 

Points in a skewed structured grid still exhibit con- 
stant spacing in x and y, however, at each level, a 
constant shift in x from the previous level results in 
a grid of parallelograms. The amount of skew in a 
skewed structured grid is 0 , the angle formed by a 
perpendicular to one of the quadrilateral sides at a 
corner (Figure 3). 

Definition 6 The Delaunay angle cut-off is the an- 
gle beyond which a skewed structured grid is no longer 
Delaunay-embeddable. 

Definition 7 A simple stretched structured grid is a 
structured grid where the location of each is given 
by: 
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Xij = i ■ 8 + ^2dx(l + e) k 
k= o 

. = r °. j = o 

J l ££i *(! + «)*. 3 = 1>2, . . n 

where dx < s/2, 0 < e < s, 0 < i < m, and 0 < j < 
n. 

The simple stretched structured grid is one which 
exhibits a constant multiplicative growth in spacing 
between levels. This will be referred to as a stretched 
structured grid for the remainder of the paper. 

The following facts are to be noted: The center of 
any circle which passes through Pi = (—s/2,0) and 
P 2 = (s/2, 0) will lie on the line x = 0 (in general, the 
center will lie on the line which is the perpendicular 
bisector to line setment PiP 2 ). Delaunay triangula- 
tions have the following properties: For any convex 
quadrilateral in the triangulation, the diagonal is se- 
lected such that the minimum angle is maximized; 
the diagonal selected is not necessarily the shortest 
diagonal. For the degenerate case of four co-circular 
points, one of two edges may be selected; in such 
cases, the edge which is a member of the edge set E 
is chosen. 
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3 Skewness and Aspect Ratio 

Given a set of quadrilaterals, ideally the principle 
direction of the skewed triangles should follow the 
original boundaries of the quadrilaterals. However, 
if there is skew (i.e., the quadrilaterals are parallelo- 
grams rather than rectangles), then it is possible for 
the Delaunay triangulation to break the quadrilateral 
boundaries. This section describes the conditions un- 
der which this happens. 

The case of a skewed structured grid with s > A is 
studied (Figure 4). The goal is to produce a Delaunay 
triangulation of the points in the vertex set V such 
that each triangle lies between the lines y = jh and 
y = (j + 1)A. In other words, we will determine the 
restrictions on the grid such that the skewed struc- 
tured grid is Delaunay-embeddable. 

First let us consider any parallelogram with aspect 
ratio s/h, with $ > A. The short diagonal creates two 
angles, a and /?, with 0 opposite the side of length s . 
As the parallelogram is skewed by 0, an amount 6 is 
added to two opposing corners and subtracted from 
the remaining two corners, causing a shift of dx to 
the upper two corner points (Figure 5). 

Lemma 1 As 0 increases, a and 0 both increase . 

Proof: Consider points P 3 and P 4 at the 
upper corners of a rectangle. Shift these 
points a distance dx relative to Pi and P 2 . 

This is the same as adding 6 to angle Pi P3P4 
and angle PiP 2 P 4 . Originally, the diagonal 
from P3 to P 2 made angles a and 0 with 
the sides PiP 2 and P 2 P 4 . Since we started 
with a rectangle, a and 0 are guaranteed 
to be less than 90 degrees. Since s > A, 
IP3P2I > IP1P3I. When the points P 3 and 
P 4 are shifted a distance dx to new points 
P3 and P 4 , two new angles a * and /?' are 
achieved. Let 0i=the angle P 3 P 2 Pg and 
0 2 =the angle P 4 P 2 P 4 . a' = a + 0i and 0' = 

0 — 8 1 -i- 0 2 . It is immediately apparent that 
a' > a. 0 f > 0 because IP 3 P 2 I > IP 1 P 3 I 
implies 0 2 > 0 X while dx < s/2. □ 

Now consider adjacent quadrilaterals PiP 2 PsP 4 
and P 3 P 1 P 5 P 6 within a skewed structured grid. A 
skewed structured grid is Delaunay-embeddable for 
the degenerate case dx = 0 since the corner points of 
any quadrilateral will be co-circular. For P 1 P 2 P 3 P 4 , 
either PiP 4 or P 2 P 3 could be chosen as diagonals. 
The following facts are to be noted: 


• dx = 0 and s > h => 0 < 90 and ZP 1 P 2 P 3 > 

• ZPiP 2 P 3 = LP4P3P2 

• As a quadrilateral PiP 2 PsP 4 is skewed, edge 
P 2 P 3 will be in the Delaunay triangulation of 
its points. 

• LP\P^P 3 is the smallest angle in the Delaunay 
triangulation of PiP 2 P 3 P 4 . 

Lemma 2 If 0 < 90 degrees, then a skewed struc- 
tured grid is Delaunay- embeddable. 

Proof: From Lemma 1, we know that 
0 increases as dx increases. It follows that 
0 approaches 90 degrees as dx increases. 

Since the lines are a constant distance h 
apart, when /?=90 degrees, the diagonal line 
P 3 P 2 becomes a perpendicular bisector for 
the line P 1 P 5 , and at this point P 3 P 2 also 
bisects the angle formed by P 1 P 2 P 5 . 

The convex quadrilateral formed by 
points P 2 PaP 5 -P 4 (Figure 6 ) has as its diag- 
onal either P3P4 or P 2 Pb- When dx = 0 , 
we know that P 3 P 4 is selected as the di- 
agonal, because ZP 4 p 3 P 2 is the larger of 
the smallest angles in the possible trian- 
gulations of this quadrilateral. As a skew 
of 6 is introduced to the grid, P5 moves 
faster than P3, so ZP3P2P5 grows faster 
than LP1P2P3 (from lemma 1 we know that 
both will increase). When /? = 90 degrees, 
the two angles are equal, and because s > A, 
ZP 4 P3P 2 < ZP5P3P4. When 0 exceeds 90 
degrees, then because ZP 3 P 2 P 5 is increasing 
faster, it becomes the larger of the smallest 
angles in the possible Delaunay triangula- 
tions of P 2 P 3 PsP 4 » and P 2 P 5 is then selected 
as the diagonal. □ 

When this occurs, the Delaunay triangulation no 
longer includes the original quadrilateral boundaries. 

Lemma 3 If s/h < 2, then a skewed structured grid 
is Delaunay- embeddable. 

Proof: At dx = 0, 0 < 90 degrees. 0 
reaches a maximum at dx = s/2. When 
s/h < 2, 0 < 90 degrees at dx = s/2, 
and the grid remains Delaunay-embeddable. 
When s/h = 2, 0 = 90 degrees at dx = s/2, 
and the grid is Delaunay-embeddable. □ 


3 


Theorem 1 As aspect ratio increases, the Delaunay 
angle cut-off for a skewed structured grid decreases. 


grid will be Delaunay-embeddable when 0 
is based on the largest aspect ratio elements 
found in the grid. □ 


Proof: From Lemma 2, we know that 
the Delaunay angle cut-off occurs when /? = 
90 degrees. At this point, the Pythagorean 
theorem gives 

s — y/s 2 — 4 h? 

dx = 2 

The Delaunay angle cut-off 9 * is defined by 



s — \fs 2 — Ah’ 1 
2h 

- 

For s/h < 2, there is no Delaunay angle 
cut-off (lemma 3). As s/h increases, 6* de- 
creases (Figure 7), □ 

Therefore, whether or not a structured grid is 
Delaunay-embeddable is dependent on aspect ratio, 
and as aspect ratio increases, the “tolerance” for skew 
decreases. As s/h gets larger, the value of the De- 
launay angle cut-off 9 * , where the Delaunay trian- 
gulation no longer includes the original quadrilateral 
boundaries, approaches 0. 

For the case of a monotonic stretched structured 
grid, the skew angle derived above is a lower bound. 

Theorem 2 As aspect ratio increases, the Delaunay 
cut-off angle of a stretched structured grid decreases. 


Proof: We need only consider two ad- 
jacent quadrilaterals at a time from the 
stretched structured grid with aspect ratios 
s/h and s/h{ 1 + e). The skew angle 6 can 
be derived from the circumcircles for the two 
interior triangles, and is defined by 


tanO 


dx_ 

h 


s — y/ s 2 — 4/i 2 — 4e/i 2 — e 2 h 2 


(2 + e)h 



■ 4 - c(4 + c) 


Since successive levels have different (in- 
creasing) values of 6 for e positive, the 


4 Conclusions 

We have shown the limitations of Delaunay triangula- 
tions of points from structured grids for aerospace ap- 
plications. For the general case of points distributed 
along fixed contours, we have shown a restriction on 
the aspect ratio for which Delaunay triangulations 
can be directly obtained. By imposing a structure 
on the point distribution, we have demonstrated the 
relationship between aspect ratio and quadrilateral 
element skew on the maintenance of contours from 
structured grids. 
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Figure 1 - Portion of a CFD Grid 
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Figure 2 - Delaunay Triangulation of CFD Grid in Figure 1 
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Figure 4 - A Skewed Structured Grid 
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